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Summary 

This report describes results of the functional flight test conducted with the Positive Systems’ ADAR 5500 
sensor system (serial number 8, linear configuration) near Winslow, Arizona on June 30 and July 1, 1999. 
The in-flight test is one component of the NASA Scientific Data Purchase (SDP) Validation and 
Verification (V&V) process. It allows to measure characteristics of the entire sensor system affected by 
both performance of the sensor during a flight and post-flight image processing. The following 
characteristics were analyzed: changes of dark digital numbers (DN’s), radiometric linearity, signal-to- 
noise ratio (SNR), spatial resolution, and geolocation accuracy. The measured characteristics were 
compared with the image product specifications defined in the Positive Systems’ SDP contract. 
Dependence of the dark DN’s on several factors was analyzed, but no significant correlation was found. 
However, the observed changes in dark DN’s were relatively small, which justifies usage of a constant 
value in the dark DN subtraction procedure during post-processing. Dependence of measured at-sensor, 
in-band radiance (in arbitrary units) on measured in-band ground reflectance is very well described by a 
linear function - The sensor fulfills the linearity requirement. Measured SNR values lower than the 
contract specifications, but accuracy of that test was possibly affected by non-uniformity of the employed 
gray-scale panels. The SNR values are generally sufficiently high for most applications. SNR can also be 
improved during standard flights by using longer exposure times. Full width at half maximum (FWHM) of 
an edge response derived line spread function was used as a measure of spatial resolution. FWHM was 
generally smaller than twice the ground sample distance (GSD), in agreement with the contract 
specifications. Accuracy of the geolocation information, which is provided for the particular images in a 
metadata file, was found to meet contract requirements as well. 

A set of images acquired during the functional flight test was also used by Positive Systems to create a 
georeferenced mosaic image. Analysis of image quality and geolocation accuracy of the mosaic were 
conducted and presented separately. The same set of images was also used to test procedures 
developed for validation of image products created after standard flights scheduled under the SDP 
program. It was found that the procedures, which are described in this report, could be applied in the 
validation process. 

The main recommendation from analysis of this functional flight test is that an exceptional care should be 
taken about size and uniformity of the gray-scale panels used in such in-flight testing. Both affect 
significantly tests of radiometric characteristics and spatial resolution. 
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Introduction 


Positive Systems, Inc., of Whitefish, Montana uses its airborne ADAR System 5500 sensors to acquire 
multispectral remote sensing images delivered to NASA under the Scientific Data Purchase (SDP) 
program. Functional flight tests are one component part of the SDP Verification & Validation process. 
Other components include laboratory and standard flight characterization. These are covered in other 
reports. The functional flight test of the Positive Systems’ ADAR 5500 sensor system with the serial 
number 8 (SN8) was conducted on June 30 and July 1, 1999 at the Hopi Reservation located southeast 
of Winslow, Arizona (see Appendix A and Appendix B). During the test, a six-step gray scale target was 
deployed on the ground at coordinates of 34° 51’ 18.865” N latitude and 110° 39’ 31.17” W longitude (see 
Figure 1) [MTL 1999], Two independent teams, one from the MTL Systems, Inc., and the other from the 
NASA Stennis Space Center (Ground Reference Information Team, GRIT), measured spectral 
reflectance of the gray panels using calibrated spectroradiometers. A network of high-precision geodetic 
targets was also located at the test site. Geolocation of the targets was measured with Global Positioning 
System (GPS) survey equipment [Shingoitewa-Honanie and Jenner, 1998] to xx precision. During the first 
day of the test, the gray-scale target was overflown several times using two mutually perpendicular 
headings and four different altitudes. Using the different altitudes, images were acquired with ground 
sample distances (GSD) of 0.3, 0.5, 0.75, and 1 meter. During the second day, images were acquired in 
seven parallel flight lines over the geodetic target range with GSD of 0.75 meter. After the flights, Positive 
Systems processed all the images as described in Appendix C and created a georeferenced mosaic from 
the images acquired on the second day. 

This documents describes results of in-flight testing of the following sensor system characteristics: 

• Dark DN’s 

• Radiometric linearity 

• Signal-to-noise ratio 

• Spatial resolution 

• Geolocation accuracy 

The set of images acquired on the second day was also used to test procedures developed for validation 
of image products acquired during standard flights. The product validation procedures and results of their 
testing are described later in this document. Evaluation of geolocation accuracy of the mosaic image is 
presented in a separate report [LMSO 1999b], 
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Figure 1. Six-step gray scale target deployed at the test site near Winslow, Arizona. 


Dark DN’s 

Dark images were aqcuired on the ground and in-flight at all the altitudes. At every altitude, two series of 
10 images each were collected. Analysis was performed separately for each series. A histogram of DN 
distribution was created for each band of every image. All the histograms in one series were used to 
calculate an average histogram and its standard deviation. Examples of the average histograms for 
different spectral bands are shown in Figure 2. 

The dark images were also used to calculate mean dark DN and its standard deviation for each series. 
Correlation of the mean dark DN’s and the average histograms with such quantities as acquisition time, 
altitude, and temperature were examined. No clear dependency on time and altitude of the acquisition 
was observed. The only significant correlation was the one with temperature (see Figure 3). Mean dark 
DN’s become smaller when temperature increases. This is inconsistent with the dark current properties of 
a Si: detector array. A possible explanation may be in thermal characteristics of the sensor electronics 
(such as amplifier offsets). Nevertheless, the mean dark DN’s changed only between 7.3 and 9.4 for the 
entire observed temperature range. Such a small change justifies usage of a constant value of 8 (instead 
of a more precise value) in the dark DN subtraction applied by Positive Systems during post-processing: 
No significant radiometric error is introduced by the approximate process. 
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c147_.5_2dark - Band 1 c147_.5_2dark - Band 2 



Figure 2. Mean histograms of dark DN distributions generated from images acquired at 
altitude suitable for 0.5 m GSD. Error bars indicate standard deviations. 


ADAR 5500 SN8 lin - in-flight dark counts 



Figure 3. Dependence of mean dark DN's on temperature. 







Radiometric linearity 

Images of the six-step gray scale target were used for in-flight characterization of the sensor radiometric 
linearity (see Figure 4). Dependence of in-band at-sensor radiance on in-band ground reflectance was 
evaluated instead of comparison of the measured radiance with the one calculated from atmospheric 
radiative transfer models. The former is more closely related to the needs of scientists who are users of 
the commercially acquired images, while the latter is prone to many uncertainties inherent to the 
atmospheric modeling [Freedman 1999], Dry atmospheric conditions of the Arizona test site may have 
contributed significantly to success of the performed analysis. 



Figure 4. Image of the six-step gray scale target acquired by 
the ADAR 5500 sensor (band 3) with GSD of 0.75 m. 

Vehicles of the ground support team are visible in the right- 
side corners of the image (compare Figure 1). 

Values of in-band radiance measured by the sensor (in arbitrary units) were extracted from the images. A 
rectangular region-of-interest (ROI) was selected for each panel area (see Figure 5), and mean and 
standard deviation of DN’s for all the pixels belonging to the ROI were calculated. The calculations were 
performed for each band separately. Mean DN was used as the measure of radiance for a given panel, 
while standard deviation was used to estimate noise (see section Signal-to-noise ratio). 

In-band reflectance of the gray panels was calculated from results of the ground measurements. While 
the measurements by the two teams produced slightly different results (see Figure 6), the spectra 
provided by MTL were used in the calculations of in-band reflectance. In-band reflectance pk is defined 
here by the following equation: 
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Figure 5. Region-of-interest selected for the radiometric 
analysis. 


P ‘ 

where X is the radiation wavelength, p(Z) is the spectral reflectance of the panel surface and R k (X) is the 
spectral response function of the k th band of the sensor. Spectral response functions measured during 
characterization of the sensor in the Commercial Instrument Validation Laboratory were used in the 
calculations [LMSO 1999a]. Both integrals were evaluated numerically with integration limits from 400 nm 
to 950 nm and with the increment of 1 nm. To have both the reflectance spectrum and the spectral 
response sampled at the same wavelengths, the functions were interpolated with splines (see Figure 7). 

Linear regression between in-band radiance and in-band reflectance was calculated using all the six gray 
panels from a single image. Gray-scale targets from different images were analyzed separately as the 
images were acquired at various times during the day under different solar illumination conditions. 
Examples of the main results of the linearity tests are shown in Figure 8. In all the cases studied, 
dependence of the radiance on the reflectance is decidedly linear as specified in the requirements. 
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Figure 7. Reflectance spectra of the gray panels and 


specirai response ot me huhr oovu ono sensor useo in 
the calculations of the in-band reflectance. 
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Figure 8. Examples of dependence of measured in-band radiance (in arbitrary units) on 
ground in-band reflectance. GSD equals 0.3 m (top left), 0.5 m (top right), 0.75 m (bottom 
left), and 1 m (bottom right). 


Sianal-to-noise ratio 

— 0 

Mean DN’s were divided by standard deviations to calculate values of signal-to-noise ratio (SNR) for the 
gray panels. Examples of dependence of the calculated SNR’s on in-band reflectance are presented in 
Figure 9. SNR monotonically increases with increasing reflectance (as well as radiance), as expected for 
a photon-limited sensor [LMSO 1999a]. For the most reflective panel, this relation is clearly violated as a 


- 12 - 






result of radiometric non-uniformity of the surface (compare Figure 1). As the entire acquisition of the 
images on June 30, 1999 lasted for a considerable length of time, solar zenith angle also changed 
significantly. The solar zenith angles were calculated using the MODTRAN program and location/time 
data provided in the GPS file. See section V of the Customer Data Sheet in Appendix C for more 
information about this file. Images, for which the six-step gray scale targets were analyzed, were acquired 
with the solar zenith angle of 13.8, 19.7, 26.4, and 35.4 degrees. These images have GSD of 0.75, 1, 0.5, 
and 0.3 m, respectively. The images with GSD of 0.5 m were collected under illumination conditions that 
are the closest to the solar zenith angle of 30° defined in the requirements. For that image, SNR is about 
60 for surface with 20% reflectance. The SNR specification for this tarp is that the SNR be greater than 
90. While this number is lower than the requirement, there is significant uncertainty in this in-flight SNR 
measurement affected by sizes of the gray panels and especially their uniformity. Moreover, the ADAR 
5500 sensor allows for selection of exposure times used in acquisition of the images. During the 
functional flight test, short exposure times were used to minimize image saturation in the gray-scale target 
area. However, during standard flights, images can be acquired using longer exposure times, which will 
increase SNR at a potential cost of more extensive saturation. 



Figure 9. Examples of dependence of signal-to-noise ratios on in-band reflectance with solar 
zenith angle of 26.4° (left) and 13.8° (right). 


Spatial resolution 

Full Width at Half Maximum (FWHM) of a line spread function (LSF) is used as a measure of spatial 
resolution of the sensor system. Line spread functions are derived from edge responses by differentiation. 
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There are three edges in the six-step gray scale target between pairs of adjacent panels. The edge 
between the most reflective panels was selected for the analysis because only that one provided 
sufficient contrast. The edge responses were measured and analyzed using a modified knife-edge 
technique [Tzannes and Mooney 1995]. A rectangular region containing a slightly tilted horizontal or 
vertical edge was extracted from an image of the gray panels as shown in Figure 10. 



Figure 10. Rectangular edge-response region selected 
from the gray panels for analysis of spatial resolution. 


In such a region, each line across the edge forms an approximate edge response. Exact edge responses 
(in the direction perpendicular to the edge) are obtained when distances are additionally scaled by cosine 
of the tilt angle. To avoid numerical differentiation, which can enhance noise, a differentiable function is 
fitted to the edge responses. In the approach presented in this report, the function has the sigmoidal form 
(known also from the Fermi distribution): 


/(*) = 


a_ 

1 + exp 


c 


b 


However, for each edge response in the selected region, position of the edge ( b ) is different because of 
the tilt. Shape of the edge causes that all the positions b are located on a straight line. Therefore, the 
fitting is performed for all the edge responses simultaneously using the formula: 


«/(*) = ' 


a 


1 + exp 


x-b x i-b 2 


+ d 


The parameters a, b u b 2 , c, and d are common for all the edge responses, while the difference in the edge 
position is introduced by the index (/). Tangent of the tilt angle is equal to the absolute value of the 
parameter b x . Figure 11 shows an example of the edge responses and the best fit to them using the 
sigmoidal functions. 
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GSD 0.5 m\c147_.5_01_271r Ian - Band 3 ( 157 167 1389 1417) GSD 0.5 m\c147_.5_01_271r.lan - Band 3 ( 157 167 1389 1417) 



Figure 11. An example of measured edge responses (left) and the best fit to them 
with sigmoidal functions (right). 


Finding the parameters b\ and b 2 is equivalent to shifting the edge responses to a single reference 
location so that all the edge points are aligned. Superimposing all the shifted edge responses creates a 
new one with a finer spatial sampling (see Figure 12). This illustrates how the modifications to the knife- 
edge method allow to overcoming the main difficulty in applying such a method for digital sensors which 
employ discrete sampling of the image space. 


B 





Figure 12. Superimposed edge responses 
and the fitted sigmoidal function. 
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The sigmoidal function with parameters obtained from the best fit is differentiated analytically to derive the 
LSF and its FWHM (see Figure 13). The LSF is sampled with a sufficient rate to avoid aliasing in the 
digital Fourier transform process, and the FFT algorithm is applied to it to compute the system modulation 
transfer function (MTF). Nyquist frequency is also marked on the MTF plot shown in Figure 14. 


GSD 0.5 m\c147_.5_01_271r.lan - Band 3 { 157 167 1389 1417) 



GSD 0.5 m\c147_.5_01„271r,lan - Band 3 ( 157 167 1389 1417) 



Figure 13. Line spread function 
derived from the edge responses 
shown in Figure 11 and Figure 12. 


Figure 14. MTF estimated from the 
LSF shown in Figure 13. 


Spatial resolution was analyzed for all the images of the gray scale target acquired during the first day of 
the functional flight test. The results for both cross-track and along-track direction are shown together in 
Figure 15. Some dependency of the results on the GSD can be noticed, but generally, FWHM is smaller 
than twice the GSD. This means that the tested sensor system, which includes both the ADAR 5500 SN8 
sensor itself and the post-processing procedures, meets the requirements for spatial resolution. 


Geolocation accuracy 

Approximate geolocation of a center point of each image is provided in the GPS file. Accuracy of those 
data was evaluated by measuring twice the distance between the image center point and another point, 
which location is known with high accuracy. One measurement was based on distance (in pixels) 
between the points and an estimated value of GSD (based on ground elevation and aircraft altitude). The 
other measurement was based on the earth's radius and the difference between geographic coordinates 
(latitude and longitude) of the points. The geodetic targets became the points with accurate geolocation. 
Seventeen images of the geodetic targets were identified. An example of such an image is shown in 
Figure 16: the target is visible as white cross on black background. Pixel coordinates of the intersection of 
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FWHM / GSD 


where: Ad- distance error 


a- instantaneous field of view of the sensor 

h - aircraft altitude 

h 0 - elevation of the geodetic target 

R - earth's radius 

Anaxj y'max - image size 

io,jo - pixel indices of the geodetic target on the image 
X, 9 - longitude and latitude of the image center 
X 0 , So - longitude and latitude of the geodetic target 



Figure 15. Ratios of FWHM to GSD determined in tests of spatial resolution based on 
measurements of edge responses. 




Figure 16. Image of a geodetic target. 


Results of the calculations of geolocation accuracy are shown in Figure 17. The sensor system meets 
requirements for the geolocation accuracy because about 90% of the points (15 of 17) have location 
errors smaller than 100 m. The mean error is 49 m with standard deviation of 37 m. 


Accuracy of GPS data 



Figure 17. Distance errors measured during evaluation 
of geolocation accuracy. 
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Tests of standard flight validation 

Verification and Validation (V&V) activities for imagery acquired under the NASA Scientific Data Purchase 
program using the Positive Systems’ ADAR System 5500 sensors consist of three parts: 

• Laboratory characterization of the sensors 

• In-flight testing of sensor performance 

• Validation of image products 

Initial tests are performed in a laboratory environment to characterize spectral, spatial and radiometric 
response of every ADAR System 5500 sensor used in the Data Purchase program. Results of such tests 
have already been reported (see, for example, LMSO 1999a). Further characterization comes from 
functional flight tests during which the sensors acquire images over well-characterized spatial, 
radiometric, and geodetic targets. Previous sections of this report describe results of such in-flight testing. 
Finally, all images produced during standard flights, scheduled under the Data Purchase program, are 
validated to ensure consistent quality of image acquisition and post-processing. Images acquired on the 
second day of the testing at the Hopi site can be considered as a model of imagery collected during a 
standard flight. Therefore, these images were used to test procedures used for validation of image 
products delivered to NASA by Positive Systems under the name IM-R1I-55, Imagery with Post-Flight 
Correction. Results of applying the validation procedure to this set of images are presented later in this 
section. The following paragraphs present the procedure as it may be described in work instructions. 

The objective is to ensure that general image quality, overlaps between images, spectral (band-to-band) 
registration, and possibly percentage of cloud cover as well as spatial resolution, conform to the contract 
requirements. Geographic coverage of a complete data set is verified during the Verification of 
Commercial Data Shipments process (see CRSP-WI-22). 

General Image Quality 

The images are evaluated to ensure proper sensor operation during the standard flight (i.e., no shutter 
failures, no significant amount of “bad” pixels, adequate settings of exposure time, etc.). The evaluation is 
based on distributions of DN’s in each band of an image. As it would be extremely difficult for product 
validation personnel to generate and analyze interactively histograms of the DN’s distributions for every 
image in a data set, the distributions are characterized by their first four moments (mean, standard 
deviation, skewness, and kurtosis) calculated using custom software (see Appendix E) and presented in a 
graphical form. Additional statistics, which are calculated for each band, include minimum DN, maximum 
DN, and percentage of pixels with the saturated DN value (the largest DN possible for a given 
quantization). Product validation personnel use the image statistics to detect anomalies in the image 
quality: Images with atypical statistics should be examined interactively using remote-sensing image 
processing software such as ERDAS Imagine or ENVI. 
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Results 

The entire set of 82 images acquired at the Hopi site in seven parallel flight lines on July 1, 1999 was 
used to test the validation procedures. The first part of image statistics is shown on graphs in Figure 18. 
For personnel conducting the image product validation, such graphs are presented interactively, with 
color coding of bands and ability to zoom in to identify particular images. The index included in Appendix 
D allows for precise identification of the images. Mean DN’s indicate changes in solar illumination during 
the time of the flight. Effects of bi-directional reflectance may also be visible. The mean values are quite 
low, especially for band 1, which indicates that exposure times were slightly too short, but that was 
required by the high priority of avoiding saturation of images of the gray scale target. Standard deviations 
of DN distributions are also affected by that. Images that include the gray scale target have significantly 
larger skewness and kurtosis. 


c147b_blk - mean c147b_blk - standard deviation 




Image # 


Image # 


c147b_blk -skewness c147b blk- kurtosis 




Figure 18. First four moments of the DN’s distributions for every image in a data 
set! mean Iton left), standard deviation (ton rinhtV skewness {bottom leftT and 

% r # » — V l- ' /» V — ft 

kurtosis (bottom right) 

Minimum and maximum DN’s are shown in Figure 19. The graph of minimum DN’s indicates that almost 
in every image there are pixels with DN equal to zero, at least in some bands. These pixels appear along 
edges of the images, in the corner areas. This is an artifact created during the spectral registration 
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processing applied by Positive Systems: the images were not sufficiently trimmed to remove the “empty” 
pixels. The graph of maximum DN’s shows that only a small number of images have saturated pixels. 


c147b bik - minimum c147b_blk - maximum 



Image # Image # 

Figure 19. Minimum and maximum DN's in each band of every image in a data set. 

Figure 20 shows that there are at most just a few saturated pixels in the images. Although the saturation 
somewhat coincides with the presence of the gray scale target in the images, close examination with the 
ENVI software showed that the saturation is caused by reflections from other objects than the gray panels 
and does not affect the radiometric tests. All the image statistics suggest that quality of images in this 
data set meets expectations. 



Image # Image # 


Figure 20. Portion of saturated pixels. Figure 21. Number of regions selected in 

each image to evaluate spectral registration. 


The number of regions selected in each image for evaluation of the spectral registration is shown in 
Figure 21. For all the images, the number of regions exceeds 400. It indicates that the program did not 
have problems with calculating the spectral registration in this data set. Mean misregistration is about 
0.35 pixel with standard deviation of about 0.2 pixel, for all the images in the data set (see Figure 22). 
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Figure 23 shows that portion of regions misregistered by more than 0.5 pixel is approximately between 
10% and 20%, while the portion of regions misregistered by more than one pixel is always smaller than 
2%. Overall, in the entire data set, less than 0.4% of the tested regions are misregistered by more than 1 
pixel. 


c147b_blk - band-to-band registration c147b_blk - band-to-band registration 




Figure 22. Mean and standard deviation of spectral misregistration. 


c147b_blk - band-to-band registration c147b_blk - band-to-band registration 



Figure 23. Portion of regions with spectral misregistration larger than 0.5 pixel f 

(left) and one pixel (right). | 
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Calculations of the image overlaps generated the results with the following statistical values: 


75 endlaps 
Mean endlap 
Std. dev. 
Minimum 
Maximum 


42.278194% 

2.464135% 

35.397453% 

45.769517% 


Portion of endlaps smaller than 35% 


0 . 000000 % 


70 sidelaps 
Mean sidelap 
Std. dev. 
Minimum 
Maximum 


43.056960% 

2.668369% 

38.276698% 

47.912062% 


Portion of sidelaps smaller than 35% 


0 . 000000 % 


All the endlaps are in the range from 35% to 46%, and all the sidelaps are in the range from 38% to 48%. 
Therefore, all the overlaps are larger than the required 35%: they perfectly meet the requirements. 


The calculated extent of the particular endlaps and sidelaps is listed below. The marks "»" or "<>" at the 
beginning of each line mean that the images were acquired in the same or in the opposite direction, 
respectively. 
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"02" 

”018 

37.961% 

» 

cl4 7b" 

blk" 

"02" 

~018 

cl47b" 

blk 

'02" 

“019 

38.2708% 

>> 

cl 4 7b" 

blk" 

"02‘ 

"019 

cl 4 7b" 

blk 

’02" 

"020 

38.0717% 

» 

cl4 7b~ 

"blk 

"02" 

“020 

cl47b" 

"blk" 

'02' 

’021 

37.8686% 

» 

Cl4 7b" 

blk" 

"02 _ 

"021 

cl 4 7b" 

"blk 

"02" 

"022 

38.1171% 

» 

cl 4 7b" 

blk" 

"03~ 

"023 

cl47b" 

"blk 

"03" 

"024 

43.1059% 

» 

cl47b" 

"blk 

"03~ 

"024 

cl 4 7b" 

"blk" 

03~ 

‘025 

42.4066% 

» 

cl47b" 

’blk" 

"03~ 

"025 

cl47b" 

"blk 

"03“ 

026 

42.1948% 

» 

cl47b 

'blk 

’03~ 

"02 6 

: cl47b 

"blk 

03“ 

"027 

42.1117% 

» 

cl47b” 

"blk 

"03~ 

"027 

: cl47b_ 

"blk 

03_ 

[028 

42.0475% 


*1 A *~I L. 

-i i ~ 


f\ 0 n 

1 A *~3 U. 

”i_ -j 

r\ n 


a -i nmo. 


UX*i IU 

u± jC 

VJ 

O 

LJX / JJ 

ui K 

uo 


*iJ L .//<£,/ 

>> 

cl47b 

"blk~ 

03~ 

‘029 

cl 47b" 

blk 

03“ 

030 

42.2414% 

» 

Cl47b 

"blk 

"03" 

030 

cl47b 

"blk 

03" 

031 

41.954% 

» 

c!47b 

‘blk 

03 

*031 

cl47b 

"blk 

03" 

032 

41.5252% 

» 

Cl47b 

blk 

03 

032 

cl47b 

"blk 

03~ 

"033 

41.5205% 

» 

cl47b 

‘blk 

03 

033 

cl47b 

"blk 

03" 

034 

40.9372% 

» 

Cl47b 

blk 

04 

035 

cl47b 

"blk 

04 

036 

45.2515% 

» 

Cl47b 

blk 

04 

036 

cl47b 

"blk 

04 

"037 

45.0603% 

» 

c!47b 

blk 

04 

037 

cl47b 

"blk 

04 

038 

45.041% 

» 

cl47b 

blk 

04 

038 

c!47b 

blk 

04 

039 

45.4522% 
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» 

cl 4 7b 

blk 

04 

039 

cl47b 

blk 

04 

040 

45.4301% 

» 

cl47b 

blk" 

“04 _ 

040 

cl47b" 

blk" 

'04' 

041 

45.1187% 

» 

c!47b 

blk" 

“04" 

~041 

cl47b" 

blk 

"04" 

042 

44.8028% 

» 

cl47b~ 

blk" 

"04" 

~042 

cl47b 

blk" 

"04* 

043 

44.2939% 

» 

cl47b 

blk 

"04" 

043 

cl47b" 

blk 

"04" 

044 

35.3975% 

» 

cl47b 

blk" 

"04 _ 

044 

cl47b 

blk 

"04~ 

045 

43.0502% 

» 

Cl47b~ 

blk" 

"04“ 

045 

cl47b 

blk] 

~04‘ 

046 

43.0196% 

» 

cl 4 7b 

blk 

"os" 

"047 

cl47b" 

blk' 

"os" 

048 

42.1373% 

» 

cl 47b" 

blk" 

"05" 

048 

cl 4 7b" 

blk" 

"os' 

049 

41.8637% 

» 

cl 47b" 

blk 

"os" 

04 9 

cl 4 7b" 

blk" 

‘05" 

050 

41.2115% 

» 

cl 4 7b" 

"blk 

"os" 

"oso 

cl 4 7b' 

blk 

'05' 

"051 

40.7031% 

» 

Cl 4 7b 

blk" 

"os" 

"osi 

cl 4 7b" 

blk' 

’05" 

052 

40.6394% 

» 

cl 47b" 

"blk 

"os" 

~052 

cl 4 7b" 

blk 

"os' 

053 

40.3825% 

» 

c 1 4 7b" 

blk" 

"os" 

“053 

cl 4 7b" 

blk 

‘os' 

'054 

39.9242% 

» 

cl 47b" 

blk 

"os" 

054 

cl 4 7b" 

blk" 

’os' 

“OSS 

39.5422% 

» 

cl 47b 

blk" 

05~ 

“OSS 

cl 4 7b" 

blk 

"os' 

056 

39.2048% 

» 

cl47b" 

blk" 

"05" 

056 

cl4 7b" 

blk 

"os' 

057 

38.4738% 

» 

cl47b" 

"blk 

"os" 

057 

cl 47b" 

blk" 

'os' 

~05 8 

37.8012% 

» 

cl47b" 

blk" 

"06 

“059 

cl 4 7b" 

blk 

0 6~ 

”060 

43.3605% 

» 

cl47b" 

blk 

"0 6" 

”060 

cl47b" 

blk' 

"oe" 

"061 

44.0976% 

» 

cl47b" 

blk 

0 6~ 

"061 

cl 4 7b" 

blk 

"06" 

“062 

45.3657% 

» 

cl47b 

blk" 

"0 e" 

"062 

cl47b" 

blk 

"0 6 

~063 

45.7695% 

» 

Cl 4 7b" 

blk" 

"oe' 

"0 63 

cl4 7b“ 

blk 

"0 e" 

~064 

45.4515% 

» 

cl 47b" 

blk" 

"06“ 

"064 

cl4 7b" 

"blk" 

‘oe' 

065 

45.2313% 

» 

cl 4 7b" 

blk" 

"06" 

"065 

cl47b" 

'blk' 

"06 

066 

45.6102% 

» 

cl 4 7b" 

blk 

"06" 

"066 

cl47b" 

blk 

"0 6 

067 

44.9615% 

» 

cl 4 7b" 

blk" 

0 6~ 

0 67 

cl47b" 

blk 

’06' 

"068 

44.0193% 

» 

cl 47b" 

blk" 

"06" 

"068 

c!4 7b" 

"blk" 

'oe' 

~069 

43.5026% 

» 

Cl 4 7b" 

blk" 

"06‘ 

069 

cl4 7b" 

"blk" 

'oe' 

070 

43.0087% 

» 

cl 4 7b] 

"blk" 

07~ 

~071 

cl47b" 

blk 

‘07“ 

~072 

42.1638% 

» 

cl4 7b" 

blk" 

"07 _ 

072 

c!47b" 

'blk 

’07' 

”073 

42.9481% 

» 

cl 4 7b" 

blk 

"07 _ 

"073 

cl4 7b" 

blk 

"07" 

074 

44.2406% 

» 

cl47b" 

blk 

07" 

“074 

cl 4 7b" 

"blk 

‘07' 

075 

44.6376% 

» 

cl47b" 

"blk 

"07~ 

075 

cl 4 7b] 

blk 

"07" 

"076 

44.8478% 

» 

cl47b" 

blk 

"07" 

076 

cl47b" 

blk 

*07' 

077 

45.1167% 

» 

cl47b] 

blk" 

"07" 

Oil 

cl 4 7b" 

"blk" 

"07" 

~078 

44.7489% 

» 

Cl47b] 

blk 

"07" 

~078 

cl4 7b" 

"blk 

'07' 

~079 

45.5187% 

» 

Cl47b" 

"blk" 

"07“ 

"07 9 

cl 4 7b' 

'blk 

'07" 

"080 

45.3268% 

» 

cl 47b" 

“blk" 

"07“ 

"080 

cl 4 7b" 

blk 

"07" 

“081 

44.9679% 

» 

cl47b~ 

_blk_ 

]07] 

]081 

cl 4 7b] 

_blk_ 

]07] 

~082 

44.3203% 

Sidelaps 









>< 

cl47b 

blk 

01 

001 

cl47b 

blk 

02 

022 

47.632% 

>< 

Cl47b" 

blk 

"oi" 

"002 

cl47b“ 

'blk 

~02‘ 

"022 

47.215% 

>< 

cl 4 7b' 

blk" 

"oi" 

"003 

cl 4 7b" 

blk 

"02“ 

"021 

46.7621% 

>< 

cl 4 7b' 

blk" 

"or 

"00 4 

cl 4 7b" 

'blk 

'02' 

”020 

46.5398% 

>< 

Cl 4 7b" 

“blk" 

'01" 

'oos 

cl4 7b" 

"blk 

"02" 

"019 

46.5384% 

>< 

Cl 4 7b" 

"blk" 

"or 

"006 

cl4 7b" 

"blk 

‘02' 

”01 8 

46.676% 

>< 

Cl 4 7b] 

bik_ 

'or 

"00 7 

cl4 7b" 

"blk 

’02" 

“017 

46.0042% 

>< 

cl 4 7b" 

blk" 

"or 

"OO 8 

cl47b" 

"blk 

02" 

016 

44.1818% 

>< 

cl 4 7b" 

"blk" 

"oi" 

"009 

cl47b" 

"blk 

"02" 

”015 

42.3052% 

>< 

cl 4 7b" 

blk 

’01" 

"OlO 

cl47b" 

"blk 

"02" 

014 

41.1496% 

>< 

cl 4 7b" 

"blk" 

'01“ 

'011 

cl47b" 

"blk 

"02" 

013 

40.8897% 

>< 

cl47b" 

"blk" 

‘02' 

"012 

cl47b“ 

'blk 

’03“ 

034 

39.0004% 

>< 

cl 4 7b" 

"blk 

‘02“ 

"013 

c!47b" 

'blk 

"03 _ 

033 

38.8713% 

>< 

cl 47b" 

'blk 

"02" 

“014 

cl47b" 

"blk 

"03“ 

032 

39.7743% 

>< 

cl 4 7b" 

"blk 

‘02“ 

'ois 

cl47b" 

'blk 

"03" 

031 

40.2214% 

>< 

cl4 7b~ 

"blk 

"02" 

"016 

cl47b“ 

'blk 

"03“ 

030 

40.7455% 

>< 

cl47b" 

"blk" 

'02“ 

'017 

cl 4 7b" 

'blk 

0 3~ 

029 

41.8939% 

>< 

cl47b" 

"blk" 

'02" 

"018 

cl47b" 

'blk 

’03" 

“028 

43.1004% 

>< 

cl47b" 

"blk 

"02“ 

"019 

cl47b" 

'blk 

‘03’ 

027 

43.0807% 

>< 

cl4 7b" 

"blk" 

"02" 

"020 

cl47b“ 

'blk 

"03~ 

026 

42.5728% 

>< 

Cl4 7b" 

"blk 

‘02' 

"021 

cl47b" 

'blk 

03 

025 

42.44% 

>< 

cl 47b" 

"blk 

'02“ 

‘022 

cl47b" 

‘blk 

’03' 

024 

43.2286% 

>< 

cl47b" 

"blk 

"03" 

‘023 

cl47b“ 

'blk 

"04“ 

"04 6 

46.398% 

>< 

cl47b" 

"blk 

‘03“ 

~02 4 

cl47b" 

"blk 

"04" 

"04 6 

46.7956% 

>< 

cl 4 7b" 

“blk 

"03“ 

"025 

cl47b" 

'blk 

"04‘ 

"04 5 

47.5893% 

>< 

cl47b" 

"blk 

"03" 

026 

cl47b“ 

'blk 

"04" 

"044 

47.9121% 

>< 

Cl47b" 

"blk" 

"03" 

027 

cl4 7b" 

'blk 

"04" 

“043 

47.1959% 

>< 

Cl4 7b" 

blk 

’03“ 

"02 8 

cl47b" 

"blk 

'04' 

~042 

45.6974% 

>< 

cl47b" 

“blk 

03~ 

"029 

cl47b" 

"blk 

"04" 

"041 

44.9477% 
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>< cl47bjblk_03_030 : cl47b_blk_04_040 : 44.7474% 
>< cl 4 7b_bl k_0 3_0 3 1 : cl47b_blk_04 _038 : 43.8045% 
X cl4 7b_blk__03_032 : cl4 7b_blk_04_037 : 44.1704% 
X cl47b_blk_03_033 : cl4 7b_blk_04 JD36 : 45.3886% 
>< cl47b_blk_03__034 : cl4 7b_blk_04_035 : 45.5823% 
X cl47b_blk_04_035 : cl47b_blk_05_058 : 42.6954% 
X c!47b_blk_04_036 : cl47b_blk_05_057 : 42.4712% 
x cl47b_blk_04_037 : cl47bj>lk_05_056 : 43.0655% 
X cl47b_blk_04_038 : c!47b_blk_05_055 : 43.637% 
X cl47b_blk_04_039 : cl47b_blk_05_055 : 43.6565% 
X cl47b_blk_04_040 : cl47b_blk_05_054 : 43.2055% 
X cl47b_blk_04_041 : cl47b_blk_05_053 : 42.372% 
>< cl47b_blk_04_042 : cl47b_blk_05_052 : 41.3148% 
X cl47b_blk_04_043 : cl47b_blk_05_051 : 40.5478% 
X cl4 7b_blk__04_04 4 : cl47b_blk_05_050 : 40.0714% 
X cl47b_blk_04_045 : c!47b_blk_05_049 : 39.3262% 
X cl47b_blk_04_046 : cl47b_blk_05_048 : 38.3691% 
>< c 1 4 7b_bl k_0 5_0 4 7 : cl47b_blk_06 __070 : 41.8428% 
X c 1 4 7b_bl k_0 5_0 4 8 : cl47b_blk_06_069 : 41.0685% 
>< cl47b_blk_05_Q4 9 : cl47b_blk_06_068 : 39.9459% 
X cl47b_blk_05_050 : cl4 7b_blk_06__0 67 : 40.1871% 
>< cl47b_blk_05_051 : cl47b_blk_06_066 : 40.3192% 
X cl47b_blk_05JD52 : cl4 7bjblk_06__065 : 40.1085% 
X cl4 7b_blk_05_053 : cl47b_blk_06_064 : 38.6513% 
>< c 1 4 7b_b 1 k_0 5_0 5 4 : cl47b__blk_06_063 : 38.2767% 
X cl47b_blk_05_055 : cl47b_blk_06_062 : 38.7231% 
>< cl47b_blk_05_056 ; cl47b_blk_06_061 : 39.9315% 
X c 1 4 7b_b 1 k_0 5_0 5 7 : cl47b_blk_06_060 : 41.1773% 
X cl4 7b_blk_05_058 : cl47b_blk_06__059 : 42.4058% 
X cl47b_blk_06_059 : cl47b_blk_07_082 : 45.3231% 
X cl47b_blk_06_060 : cl47b_blk_07_081 : 44.9186% 
X cl47b_blk_06_061 : cl47b_blk_07_080 : 45.2521% 
X cl47b_blk_06_062 : cl47b_blk_07_079 : 46.039% 
X cl47b_blk_06__063 : cl47b__blk_07_078 : 46.0979% 
X cl 4 7b_blk_0 6_0 64 : cl47b__blk_07_077 : 45.024 6% 
>< cl4 7b_blk_06_065 : cl4 7b_blk_07__07 6 : 43.502% 
>< c 1 4 7b_bl k_0 6_0 6 6 : cl4 7b_blk_07_075 : 42.3547% 
X cl47b_blk_06_067 : cl47b_blk_07_074 : 41.7012% 
>< cl47b_blk_06_068 : cl47b_blk_07_073 : 42.2176% 
>< c 1 4 7b__blk_0 6_0 6 9 : cl 4 7b_blk_07__072 : 43.1648% 
X c!47b blk 06 070 : cl47b blk 07 071 : 43.9692% 
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Appendix A 


FLIGHT LOG 
(6-30-99) 


Customer name: 
Customer number: 
Task Order: 

Task Request: 
Flight Crew: 

Data set: 

Site Name: 


NASA - SDB 
cl47 

Task Order #1 

Task Request - Test Flight 
Willie & Brad 
a 

Winslow, AZ 


- ADAR System: 


SN-8 




- Camera 

orientation: 


1 & 2 

Inverted 



- Flight 

Altitude: 


9,100, 

11,700, 

14,900, 18 

!, 100 'MSL 

- Ground 

elevation: 


5,300' 




- Wheels 

up at: 


11:20 

MDT 



- Weather 

forecast : 


Sunny 




- Sky conditions on site: Sunny 




- On station at: 


12:00 

MDT 



EXPOSURES 

; 






Band 1 

ISO=200, 

f 

2.8, shutter speed= 

1250 


Band 2 

ISO=200, 

f 

2.8, shutter speed= 

2500 


Band 3 

ISO=200, 

f 

4.0, shutter speed- 

800 


Band 4 

ISO=200, 

f 

5 . 6, shutter speed= 

800 


NOTE: Data is very dark in bands 1 & 

[ 2 due tc 

i the fact 

that we could not 

any saturation of the 

radiometric panels. 



LINE# 

HEADING 


IMAGE# 1 s 

NOTES 



grnd_dark 

n/ a 


001-010 

Dark imgs 

on the ground. Temp 34. 8C 

. 75_ldark 

n/a 


013-023 

1st dark 

at .75mpp. 

Temp 22. 7C. 

Single Lines: .75mpp 






.75 01 

083 


122-131 




.75 02 

353 


132-142 




. 75 03 

083 


143-152 




.75 04 

353 


153-163 




.75 05 

083 


164-172 




.75 06 

353 


173-183 




.75 2dark 

n/a 


184-193 

2nd dark 

at .75mpp. 

Temp 16. 9C. 

Single Lines: Impp 






1 ldark 

n/a 


194-203 

1st dark 

at Impp . 

Temp 16.4C. 

1 01 

353 


204-210 




1 02 

083 


211-217 




1_03 

353 


218-224 




1 04 

083 


225-231 




1 _ 05 

353 


232-238 




1_06 

083 


239-245 




1 2dark 

n/a 


246-255 

2nd dark 

at Impp. 

Temp 19.3C. 
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Single Lines: . 5mpp 


.5 ldark 

n/a 

256-265 

1st dark at 

. 5mpp 

Temp 

20. 7C 

.5 01 

083 

266-276 





.5 02 

353 

277-289 





.5 03 

083 

290-301 





.5 04 

353 

302-314 





.5 05 

083 

315-326 





. 5_0 6 

353 

327-338 





.5 2dark 

n/a 

339-348 

2nd dark at 

. 5mpp 

. Temp 

25. 0C 

Single Lines 

; : . 3mpp 






.3 ldark 

n/a 

349-369 

1st dark at 

. 3mpp 

Temp 

27. 2C 

.3 01 

353 

370-382 

rough flight 

at 0 

. 3mpp . 


.3 02 

083 

383-393 





.3 03 

353 

413-420 

reflight of 

line i 

03. 


.3 04 

083 

421-428 





.3 05 

353 

429-439 





• 3_06 

083 

440-447 





.3 2dark 

n/a 

448-455 

2nd dark at 

. 3mpp 

. Temp 

31. 1C 

Last image: 


15:30 MDT 
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Appendix B 


FLIGHT LOG 
(7-1-99) 


Customer name: 

Customer number: 

Task Order: 

Task Request: 

Flight Crew: 

Data set: 

Site Name: 

ADAR System: 

Camera orientation: 
Flight Altitude: 

Ground elevation: 

Wheels up at: 

Weather forecast: 

Sky conditions on site: 
On station at: 


NASA 

cl47 

Task Order #1 

Task Request - Test Flight 
Willie & Brad 
b 

Winslow, AZ 
SN-8 

1, 2 inverted 
14,900ft MSL 
5, 300ft 
10:30 MDT 
Sunny 
Sunny 
11:00 MDT 


EXPOSURES: 




Band 1 

ISO=200, f 

2.8, 

shutter speed= 1250 

Band 2 

ISO=200, f 

2.8, 

shutter speed= 2500 

Band 3 

ISO=200, f 

4.0, 

shutter speed= 800 

Band 4 

ISO=200, f 

5.6, 

shutter speed= 800 

LINE# 

HEADING 

IMAGE #'s N< 

blk 01 

090 


01-11 

blk 02 

270 


12-22 

blk 03 

090 


23-34 

blk 04 

270 


35-46 

blk 05 

090 


47-58 

blk 06 

270 


59-70 

blk_07 

090 


71-82 

Last Image: 


11:30 

MDT 


NOTES 
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Appendix C 


Customer Data Sheet 

RE: c147 - Delivery Task Order #1 ; Task Request - Test Flight; 

Site: Winslow, A Z 

I. Image Format 

A. Images are in ERDAS “.Ian” format 

B. This delivery includes images captured in two datasets: 


c147 

flight day 1 

c147b 

flight day 2 

naming convention: 

c147*_xx_ 

_??_###r.lan 

c147 

project ID 

* 

dataset distinction (b for block site ) 

XX 

imagery descriptor (.3, .5, .75, 1 , dark, blk) 

?? 

flight line 

### 

scene number 

r 

indicates that the image has been spectrally co-registered 
(not present for dark images) 

Jan 

ERDAS .Ian extension 



Dataset 

Images 

# of Images 


flight day 1 

0.3 mpp GSD 

59 



0.5 mpp GSD 

73 



0.75 mpp GSD 

62 



1 mpp GSD 

42 



dark images 

103 


flight day 2 (b) 

block site 

82 

D. 

Number of images delivered: 

421 


E. 

Data volume: 

2.52 GB 


Flight Line Layout 



A. 

Average Spatial Resolutions: 

0.3, 0.5, 0.75, and 1 meter/pixel GSD 

B. 

Desired Overlap: 

35% end and 35% side for block site 

C. 

Aircraft Altitude: 

0.3 mpp 

9,100’ MSL 



0.5 mpp 

11,700’ MSL 



0.75 mpp 

14,900’ MSL 



1 mpp 

18,100’ MSL 



block site 

14,900’ MSL 



dark images 

acquired at 9,100’, 11,700’, 
14,y00\ and 18,100 Msi. 

D. 

Average Image Footprint Size: 

300 x 450 meters @ 0.3 meter/pixel GSD 
500 x 750 meters @ 0.5 meter/pixel GSD 
750 x 1 125 meters @ 0.75 meter/pixel GSD 
1000 x 1500 meters @ 1 meter/pixel GSD 
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E. Desired Heading: 

1 . Individual tarp flights at resolutions of 0.3mpp, 0.5mpp, 0.75mpp, 
and Impp GSD: 



3. Dark images: heading varies - not applicable. 

F. Number of Flight Lines 

1. Tarp flights: 

2. Block site: 

3. Dark images: 

G. Flight log information: 

For flight log information, please refer to the individual flight logs included with this 
delivery. There are two flight logs, provided in both hardcopy and electronic format, one 
for each dataset (each day of image acquisition). The flight log names are: 

c147_flt_log.txt flight day 1 

c147b_flt_log.txt flight day 2 

III. Sensor Information - Sensor used for this project was SN8 


6 flight lines at each resolution 

7 flight lines 

2 flight lines at each resolution 



A. 


Spectral Bands and Band Widths 


Bandwidth (nmt Color 

450-515 

525-605 

630-690 

750-900 


Blue 
Green 
Red 
Near IR 


B. 


Location of sensors in the housing: 
(direction of flight towards top of page) 



The two forward sensors, in this case #1 and #2, record images which are inverted 
compared to images captured by the other two sensors, #3 and #4. 
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IV. Dark Image Information 


Image Name 
Prefixes 

Resolution 

06/30/99 

Approximate Time 
(GMT) 

Temperature (°C) 

c147_.3_1dark 

0.3 mpp 

21:48 

27.2 

cl 47 .3 2dark 

0.3 mpp 

22:27 

31.1 

c147_.5_1dark 

0.5 mpp 

21:10 

20.7 

c147 .5 2dark 

0.5 mpp 

21:42 

25.0 

c147 .75 Idark 

0.75 mpp 

18:47 

22.7 

cl 47 .75 2dark 

0.75 mpp 

20:17 

16.9 

cl 47 1 Idark 

1 mpp 

20:28 

16.4 

c147_1_2dark j 

1 mpp 

20:58 

19.3 

cl 47 grnd_dark 

N/A 

18:05 

34.8 


V. GPS Information 


The GPS files are in electronic format on the enclosed CD-ROM and on the 8mm image delivery 
tape: 

c147_.3_gps.txt 

c147_.5_gps.txt 

c147_.75_gps.txt 

c147_1_gps.txt 

cl 4 7_dark_gps. txt 
cl 4 7b_block_gps. txt 

The GPS text file includes one world geodetic coordinate for each image captured. These 
latitude/longitude coordinates are generated utilizing the NAD 27 datum to reflect the approximate 
camera exposure station for each image. These files include scene number, dataset, direction of 
flight, exact time of image capture, approximate latitude/longitude, and elevation of the aircraft at 
the time of image capture. See the enclosed GPS technical note for a description of the GPS text 
data, provided in both hardcopy and electronic format. The electronic file on the CD-ROM is 
named gps_tech.txt. 

VI. Post Processing 


A. 


B. 


Fully Processed Images (individual tarp images and block site): 

1 . All fully processed images are on the enclosed 8mm tape. 

2. All fully processed images are in ERDAS ".Ian” format. 

3. Bands 1 & 2 have been rotated 1 80 degrees to match bands 3 & 4. 

4. Images have been corrected for sensor nonuniformities. 

5. Images have been corrected for lens vignetting effects. 

6. Images have had a dark sample subtraction value of 8 applied. 

7. Images have been spectrally co-registered using the nearest neighbor algorithm. 


Dark Images: 

1 . All dark images are on the enclosed 8mm tape. 

2. All dark images are in ERDAS “.Ian” format. 

3. For a!! dark images, bands 1 & 2 have not been rotated ISO degrees to match bands 
3 & 4. 


4. No post-processing has been performed on the dark images (i.e., no correction for 
sensor nonuniformities, no correction for lens vignetting effects, no dark sample 
subtraction, and no spectral co-registration). 
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Appendix D 

Image file index used on the horizontal axes of the graphs created during validation of image products: 


1: cl47b_blk_01_001r . Ian 
2 : cl47b m blk_01_002r . Ian 
3: cl47b_blk_01_003r . Ian 
4: cl47b_blk_01_004r . Ian 
5: cl47b_blk_01_005r . Ian 
6: cl47b_blk_01_006r . Ian 
7: cl4 7b_blk_01_007r . Ian 
8: cl47b_blk_01_008r.lan 
9: cl47b_blk_01_009r . Ian 
10: cl47b_blk_01_010r . Ian 
11: cl47b_blk_01_011r . Ian 
12: cl47b_blk_02_012r . Ian 
13 : cl4 7b_blk_02__013r . Ian 
14: cl47b_blk_02_014r . Ian 
15: cl47b_blk_02_015r . Ian 
16: cl47b_blk_02_016r . Ian 
17: cl4 7b_blk_02_017r . Ian 
18: cl47b_blk_02_018r . Ian 
19: c!47b_blk_02_01 9r . Ian 
20: cl47b_blk_02_020r . Ian 
21: cl47b_blk_02_021r . Ian 
22: cl47b_blk_02_022r . Ian 
23: cl47bjolk_03_023r.lan 
24: cl47b_blk__03_02 4r . Ian 
25: cl47b_blk_03JD25r.lan 
26: cl47b_blk_03_026r.lan 
27: cl47b_blk_03_027r . Ian 
28: cl4 7b_blk_03__028r . Ian 
29: cl47b_blkj)3_029r.lan 
30: cl4 7b_blk__03__030r . Ian 
31: cl47b_blk_03_031r . Ian 
32: cl47b__blk_03_032r.lan 
33: cl47b_blk_03_033r . Ian 
34: cl47b_blk_03_034r . Ian 
35: cl47b_blk_04_035r . Ian 
36: cl47bjolkj)4_036r.lan 
37: cl47b_blk__04__037r . Ian 
38: cl47b_blk_04_038r . Ian 
39: cl47b_blk_04_039r . Ian 
40: cl47b_blk_04_040r.lan 
41: c!47b blk 04 041r.lan 


42 : cl47b_blk_04_042r . Ian 
43: cl47b_blk_04_043r . Ian 
44: cl47b_blk_04_044r . Ian 
45: cl47b_blk_04_045r . Ian 
46: cl4 7b_blk_04_04 6r . Ian 
47: cl47b_blk_05_047r . Ian 
48: cl47b_blk_05_048r.lan 
49: cl47b__blk__05_04 9r . Ian 
50: cl47b_blk_05_050r.lan 
51: cl47b_blk_05_051r . Ian 
52: cl47b_blk_05_052r.lan 
53: cl47b_blk_05_053r . Ian 
54: cl47b_blk_05_054r.lan 
55: cl47b_blk_05_055r . Ian 
56: cl47b_blk_05_056r . Ian 
57: cl47b_blk_05_057r . Ian 
58: cl47b_blk_05_058r . Ian 
59: cl47b_blk_06_059r.lan 
60: cl47b_blk_06_060r . Ian 
61: cl47b_blk_06_061r . Ian 
62 : cl4 7b_blk__06__062r . Ian 
63: cl47b_blk_06_063r . Ian 
64: cl47b_blk_06_064r . Ian 
65: cl4 7b__blk__0 6__0 65r . Ian 
66: cl47b_blk_0 6_066r . Ian 
67: cl4 7b__blk_06_0 67r . Ian 
68 : cl47b__blk__06_068r . Ian 
69: cl47b_blk_06_069r . Ian 
70: cl47b_blk_06_070r . Ian 
71: cl47b_blk_07_071r . Ian 
72: cl 4 7b_blk_07__072r . Ian 
73: cl47b_blk_07_073r . Ian 
74: cl47b_blk_07_074r . Ian 
75: cl47b_blk_07_075r . Ian 
76: cl47b_blk_07_076r . Ian 
77: cl47b_blk_07_077r . Ian 
78: cl47b_blk_07_078r . Ian 
79: cl47b_blk_07_079r . Ian 
80: cl47b_blk_07_080r . Ian 
81: cl47b_blk_07_081r . Ian 
82: cl47b blk 07 082r.lan 



f 


Appendix E 

MATLAB function that calculates image statistics and spectral registration for a set of images: 

function [] = calc_PosSys ( dirname ) 

% Performs calculations for validation of a set of Positive Systems' ADAR 5500 imagery 

% Set data set identification name by extracting the current subdirectory name from the full 
directory name 

delimiter = fullfile( ' l f , '3* ); delimiter = delimiter (2) ; % find separator 

k = findstr{ dirname, delimiter ); k = max ( k ); % find position of the last separator 

ident = dirname (k + l:end); % extract the last subdirectory 

% Select image file format (set for ADAR 5500 imagery) 

filetype = 'LAN'; 

file_ext - lower ( filetype ); 

% Get list of files 


filelist = dir( fullfile ( dirname, [ ’ * . ' file_ext ] ) ); 
nframes - size( filelist, 1 ); 

% Print the image file names into the index file 

filename = fullfile ( dirname, [ 'index_' ident ' . txt * ] ); 
fid - fopen( filename, 'wt' ); 
for id = 1 : nframes 

fprintf ( fid, f %d : %s\n', id, filelist (id) .name ); 
end 

f close ( fid ); 

% Get image specifications from the first file 

filename = fullfile ( dirname, filelist (1) .name ); 
image_cube - open_image_cube_file ( filename, filetype ); 
close_image_cube_file ( image_cube ) 

% Set parameter for calculations of statistics 


sat_level = 2 A ( 8 * image_cube . nbytes ) - 1; % maximum DN (saturation level) 

% Set parameters for calculations of spectral registration 

n * 17; % cross-correlation block size 
m = 8; % peak block size 

k = 16; % resolution enhancement factor 

overlaps =0.35; % frame overlaps 

skip = 5; % no. of control regions to skip 

% Initialize matrices for statistics 


avg^frame 
std_frame 
skw_ frame 
kur_frame 
min__frame 
max_frame 
sat frame 


zeros 

zeros 

zeros 

zeros 

zeros 

zeros 

zeros 


nframes, 

nframes, 

nframes, 

nframes, 

nframes, 

nframes, 

nframes. 


image_cube . nbands 
image__cube . nbands 
image_cube . nbands 
image_cube . nbands 
image_cube . nbands 
image_cube . nbands 
image_cube . nbands 


) ; 
) ; 
) ; 
) ; 
) ; 
) ; 
) ; 


% Initialize matrices for spectral registration 


avg_misreg = zeros ( 
std_misreg = zeros ( 
len_misreg = zeros ( 
num_misreg = zeros ( 
one_misreg = zeros ( 


nframes, 

nframes, 

nframes, 

nframes, 

nframes, 


1 ) 
1 ) 
1 ) 
1 > 
1 ) 
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% Initialize matrix for band correlation 

corr_band = zeros ( nframes, image_cube . nbands * ( image__cube . nbands - 1 ) / 2 ); 

% Calculate statistics for each image 
for iframe = 1 : nframes 

% read image file from disk 

filename = fullfile ( dirname, filelist (iframe) .name ); 
image_cube = open_image_cube_f ile ( filename, file type ) ; 
image_cube . data = read_image_cube__tile ( image_cube, [] ); 
close_image_cube_file ( image_cube ) 

% calculate statistics for each band 
for iband = 1 : image_cube . nbands 

frame = double ( image__cube . data ( : , : , iband) ) ; 
frame = frame ( : ) ; 

avg_frame (if rame, iband) = mean ( frame ); 
std_frame (iframe, iband) = std( frame ); 
skw_f rame (if rame, iband) - skewness ( frame ); 
kur_f rame (if rame, iband) = kurtosis ( frame ); 
min_frame (if rame, iband) = min ( frame ); 
max_f rame (if rame, iband) ~ max ( frame ); 

sat_f rame (if rame, iband) = sum( frame == sat_level ) / length ( frame ); 
end 

% calculate misregistration between all the bands 

nbrows = floor ( ( 1 - overlaps ) * image__cube .nrows / ( skip * n ) ); 

nbcols = floor ( ( 1 - overlaps ) * image_cube . ncols / ( skip * n ) ); 

rowO ~ f loor ( ( image_cube . nrows - skip * n * nbrows ) / 2 ); 

colO - floor ( ( image_cube . ncols - skip * n * nbcols ) / 2 ); 

misreg = zeros ( nbrows, nbcols, image_cube .nbands * ( image__cube . nbands - 1 ) / 2 ); 

frame = [ ] ; 

for ib = 1 : nbrows 

rows = [ 1 : n ] + rowO + { ib - 1 ) * skip * n; 
for jb = 1 : nbcols 

cols = [ 1 : n ] + colO + ( jb - 1 ) * skip * n; 
for iband = 1 : image_cube . nbands 

f rame (:,:, iband) = double ( image__cube .data (rows, cols, iband) ); 
frame ( : , : , iband) « frame ( : , : , iband) - mean2 ( frame ( : , : , iband) ) ; 
end 

ipair =0/ 

for iband - 1 : image_cube . nbands - 1 

for jband = iband + 1 : image_cube .nbands 

xc = xcorr2 ( f rame (:,:, iband) , frame (:,:, jband) ); 
xc = xc(l : end - 1,1 : end - 1); 

xc = xc ( ( end - m ) / 2 + 1 : ( end + m ) / 2, ( end - m ) / 2 + 1 : ( end + m ) / 2) / 
% enhance spatial resolution 
u = fftshift ( fft2 ( xc ) ) ; 

v = cat ( 1, zeros ( ( k - 1 ) * size ( u, 1 ) / 2, size( u, 2 ) ) , u, ... 

zeros ( ( k - 1 ) * size ( u, 1 ) / 2, size{ u, 2 ) ) ); 

v - cat ( 2, zeros ( size ( v, 1 ), ( k - 1 ) * size ( u, 2 ) / 2 ) , v, ... 

zeros ( size( v, 1 ) , ( k - 1 ) * size ( u, 2 ) / 2 } ); 

yc = abs ( iff t2 ( v ) ) ; 

% estimate misregistration 
[ i, j ] = find ( yc == max ( yc ( : ) ) ); 

dist = sqrt ( (k*m/2+l- i { 1 ) ) A 2+ (k*m/2+l- j (1) ) A 2 ) / k; 
ipair = ipair + 1/ 
misreg (ib, jb, ipair) = dist; 
end 
end 
end 
end 

misreg = misreg (find ( misreg < 2 )); 
avg__misreg (if rame) = mean ( misregO) ); 
std_misreg (iframe) = std( misreg(:) ); 
len_misreg (iframe) = length ( misreg(:) ); 

num_misreg (iframe) = sum( misreg(:) > 0.5 ) / len_misreg (if rame) ; 
one_misreg (iframe) = sum( misreg(:) > 1 ) / len_misreg (if rame) ; 
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% calculate correlation between bands 
ipair = 0; 

for iband = 1 : image_cube . nbands - 1 

for jband = iband + 1 : image__cube . nbands 
ipair = ipair + 1; 

cor r_band (if rame, ipair) = corr2 ( image_cube .data iband) , image_cube . data ( : 
end 
end 

end 

% Save results 

filename ~ fullfile( dirname, [ 'matrix_' ident '.mat 1 J ); 
save ( filename, 'dirname', ... 

, avg_frame', ' std_frame', ' skw_f rame ' , ’kur_frame*, 'min_frame', 'max_frame', 

' avg_misreg ' , * std_misreg' , ' num_misreg' , ' len_misreg ' , * one_misreg ' , ... 

'corr_band' ) 


return 


jband) ) 


'sat frame' 
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Appendix F 

MATLAB function that calculates image overlaps for a set of ADAR 5500 images. 

function [] = frame_overlaps ( dirname, gpsfile, demdir ) 

% Uses a Postitive Systems GPS output file and scene information 
% to extract georeference data (location, time, heading, image size) 

% and calculate endlaps and sidelaps between images 

% Set data set identification name by extracting the current subdirectory name from the full 
directory name 

delimiter = fullfile ( ' 1 ' , '3' ); delimiter = delimiter (2) ; % find separator 

k = findstr( dirname, delimiter ); k = max ( k ); % find position of the last separator 

ident = dirname (k + l:end); " extract the last subdirectory 

% Select image file format (set for ADAR 5500 imagery) 

filetype = 'LAN'; 

file_ext = lower ( filetype ); 

% Set some constants (these would be passed in if a GUI were used) 

ft2km = 12 * 25.4 * le-6; % km/ft - in/ft * mm/in * km/mm 

radius = 6378.137; % earth radius [km] 

end_spec - 0.35; % endlap requirement 

side_spec = 0.35; % sidelap requirement 

ifov_x - 0 . 257e-3; % ADAR 5500 cross-track IFOV (sampling interval) [rad] 
ifovJL = 0 . 257e-3; % ADAR 5500 along-track IFOV (sampling interval) [rad] 

% Read in the GPS file and get the Dataset #, 

% Latitude, Longitude, Altitude, and Heading. 

iframe - 0; 
nlines = 0; 
dataset = []; 
f light_line = [ ] ; 

fid = fopen ( fullfile ( dirname, gpsfile ) , 'r' ); 

while ~feof( fid ) 
line = fgetl ( fid ) ; 

if i s empty ( line ); break; end; % there may be two empty lines at the end of the file 
iframe = iframe + 1; 

dataset - strvcat ( dataset, sscanf ( line, ' %*s %s', 1 ) ); 

% detect a new flight line and note its name to the list of the flight lines 
% (assuming that all the images from one flight line are listed together in the GPS file) 
if ( iframe == 1 ) I ~strcmp ( dataset (iframe, :) , dataset (iframe - 1,:) ) 
nlines = nlines + 1; 

flight_line = strvcat ( flight_line, dataset (iframe, : ) ); 
end 

% get scene (image) number 

line = fgetl ( fid ); 

scene_num (if rame) = sscanf ( line, ' %*s %*s %d', 1 ); 

% check order of images in the same flight line 

if ( iframe > 1 ) & strcmp ( dataset (iframe, :) , dataset (iframe - 1,:) ) & ... 

( scene_num( iframe) <= scene_num (iframe - 1) ) 

disp ( [ 'Warning: image ' deblank( dataset (iframe, : ) ) sprintf ( '_%03u', scene_num (if rame) ) 

* follows ' deblank { dataset (if rame - 1,:) ) sprintf ( '_%03u', scene_num (iframe - 1) 

) ] ) 

end 

line = fgetl ( fid ) ; 

latitude (if rame) = sscanf ( line, '%*s %f, 1 ); 
hemisphere = sscanf ( line, *%*s %*f %*s %s', 1 ); 
if strcmp { hemisphere, 'S' ) 

latitude (iframe) = - latitude (if rame) ; 
end 

longitude (iframe) = sscanf ( line, '%*s %*f %*s %*s %*s %f', 1); 
direction = sscanf ( line, '%*s %*f %*s %*s %*s %*f %*s %s', 1); 
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if strcmp( direction, ' W' ) 

longitude (if rame) = - longitude (if rame) ; 
end 

line = fgetl { fid ) ; 

altitude (iframe) = sscanf ( line, '%*s %f', 1); 
line = fgetl ( fid ) ; 

heading (iframe) = sscanf ( line, '%*s %f ? , 1); 
line = fgetl ( fid ) ; 
line = fgetl ( fid ) ; 

% get the frame size from the image file header 

filename = [ deblank( dataset (iframe, : ) ) sprintf ( f _%03u', scene_num (iframe) ) 'r.' file_ext 

] ; 

image_cube - read_LAN_header ( fullfile( dirname, filename ) ); 
scene_rows (iframe) - image_cube . nrows; 
scene_cols (iframe) = image_cube . ncols; 

% find ground elevation from the DEM data 

[ filename, quadrangle ] = usgsdems ( latitude (if rame) , longitude (if rame) ); 
filename = char ( filename{l} ); % convert cell to string 
[ dem, leg ] = usgsdem( fullfile( demdir, filename ), 1, ... 

[ latitude (iframe) latitude (iframe) ], [ longitude (if rame) 

longitude (if rame) ] ); 

elevation (iframe) = dem(l,l) * le-3; % [m] -> [km] 

end 

f close ( fid ) ; 
nframes - iframe; 

% Define adjacent flight lines 

adjacent_flight_line - []; 

filename = ’adjacent_flight_lines.txt'; 

if exist ( fullfile ( dirname, filename ) ) == 2 

% overwrite default line order using data in the file 
fid = fopen( fullfile ( dirname, filename ) , * r* ); 

while ~feof ( fid ) 
line = fgetl ( fid ) ; 

adjacent_flight__line = strvcat ( adjacent__flight_line, line ); 
end 

fclose ( fid ) ; 
else 

% create the default list: the adjacent flight line is the next one 

for i - 1 : nlines - 1 

adjacent_flight_line = strvcat ( adjacent_flight_line, [ flight__line (i, : ) ' ' flight_line (i + 

1 ,:) ] ); 

end 

end 

for i = 1 : nframes 

j = strmatch( [ deblank( dataset (i,:) ) ' ' ], adjacent_f light_line ); 

if isempty( j ) 
line_id(i) = 0; 
else 

line_id(i) = j; 
end 
end 

% Get half of image size [kml 

alongtrack = tan ( ifov_l * scene_rows / 2 ) . * ( altitude * ft2km - elevation ); 
crosstrack = tan ( ifov_x * scene_cols / 2 ) . * ( altitude * ft2km - elevation ); 

% Change azimuth to trigonometric angle 

angle = 90 - heading + 360 * { heading >= 270 ) ; 

% Change degrees to radians 

latitude - latitude * pi / 180; 
longitude - longitude * pi / 180; 
angle = angle * pi / 180; 

% Calculate image versors 
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vx = cos ( angle ) ; 
vy = sin( angle ); 

% Open the output file 

filename = fullfile ( dirname, [ ' endlaps_sidelaps___' ident '.txt' ] ); 
fid = fopen( filename, 'wt' ); 

% Calculate endlaps (within each flight line) 

fprintf ( fid, ' Endlaps\n' ); 


count =0; 

for i = 1 : nframes - 1 


j = i + 1; 

% ensures that only frames in the same flight line are compared 

if strcmp( dataset (i, :) , dataset (j,:) ) 
count = count + 1; 

% project distance and sizes onto the direction of mean azimuth 
if ( vx (i) * vx { j ) + vy (i) * vy(j) ) > 0 
% flight lines in the same direction 
endlap (count) = 1 - radius * ... 

abs ( { vx(i) + vx(j) ) * ( longitude (j) - longitude (i) ) * ... 

cos ( { latitude (i) + latitude (j) ) / 2 ) + 
( vy(i) + vy(j) ) * ( latitude (j) - latitude (i) ))/... 
abs ( ( vx(i) + vx(j) ) * ( alongtrack (i) * v-x(i) + alongtrack ( j ) * 


vx ( j) ) + ... 

vy ( j ) ) ) ; 

line = 


else 


{ vy (i) + vy ( j ) ) * ( alongtrack (i) * vy(i) + alongtrack (j ) 

'» ’ deblank ( dataset (i,:) ) sprintf ( *_%03u', scene_num(i) 
deblank ( dataset (j,:) ) sprintf ( '_%03u', scene_num{j) 
num2str ( 100 * endlap (count) ) '%' ]; 


% flight, lines in the opposite directions 

endlap (count) = 1 - radius * ... 

abs ( ( vx(i) - vx(j) ) * ( longitude (j) - longitude (i) ) * ... 

cos ( ( latitude (i) + latitude (j) ) / 2 ) + 
( vy(i) - vy(j) ) * ( latitude (j) - latitude (i) ))/... 

abs ( ( vx(i) - vx(j) ) * ( alongtrack (i) * vx(i) - alongtrack (j ) * 


vx(j) ) + .. 

vy ( j ) ) ); 

line = 


( vy(i) - vy(j) ) * ( alongtrack (i) * vy(i) - alongtrack (j ) * 

»>< * deblank ( dataset (i, : ) ) sprintf ( '_%03u', scene_num(i) )’:'... 

deblank ( dataset (j,:) ) sprintf { , _%03u', scene_num(j) )':'... 
num2str ( 100 * endlap (count) ) '%' ]; 


end 

if endlap (count) > 0 

fprintf ( fid, '%s\n’, line ); 
else 

fprintf ( fid, '%s !!! GAP !!!\n', line ); 
end 
end 
end 


% Calculate sidelaps (between flight lines) 
fprintf ( fid, ' \nSidelaps\n ' ); 

count = 0/ 

for i = 1 : nframes - 1 

% find the closest image from the adjacent flight line 

aist_min = Inf; 

for j = i + 1 : nframes 

% not the last line & not the same line & it is the adjacent line 

if ( line_id (i) > 0 ) & -strcmp( dataset (i, :) , dataset (j,:) ) & ... 

~isempty ( f indstr ( deblank ( dataset (j,:) ), adjacent_flight__line (line_id (i) , : ) ) ) 

% calculate distance between image centers 

dx = ( longitude (j) - longitude (i) ) * cos ( { latitude (i) + latitude (j) ) / 2 ); 

dy = ( latitude (j) - latitude (i) ); 
dd = radius * sqrt ( dx * dx + dy * dy ) ; 
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if dd < dist_min 
dist_min = dd; 
j_min = j ; 
end 
end 
end 

% if the image exists 
if dist__min < Inf 
j = j_min; 

if (vx(i) * vx(j) + vy(i) * vy(j) ) >0 
% flight lines in the same direction 
% estimate apparent endlap 
endlap_ = 1 - radius * ... 

abs ( ( vx(i) + vx(j) ) * { longitude (j) - longitude (i) ) * ... 

cos ( ( latitude (i) + latitude (j) ) / 2 ) + ... 

( vy{i) + vy(j) ) * ( latitude (j) - latitude (i) ) ) / ... 

abs ( ( vx(i) + vx(j) ) * ( alongtrack (i) * vx(i) + alongtrack ( j ) * vx(j) ) + 

( vy{i) + vy ( j ) ) * ( alongtrack (i) * vy(i) + alongtrack (j ) * vy(j) ) ); 

if endlap_ > 0 

% include sidelap only when the endlap exists 

count = count + 1; 

sidelap (count) = 1 - radius * ... 

abs ( ( vy(i) + vy(j) ) * ( longitude (j) - longitude (i) ) * ... 

cos ( { latitude (i) + latitude (j) ) / 2 ) — 


( vx(i) + vx(j) ) * ( latitude (j) - latitude (i) ))/... 

abs { ( vx(i) + vx(j) ) * ( crosstrack (i) * vx(i) + crosstrack (j ) 

vx ( j ) ) + ... 

( vy (i) + vy { j ) ) * ( crosstrack (i) * vy(i) + crosstrack (j ) 

vy ( j ) ) ) ; 

line = [ f >> ' deblank ( dataset (i,:) ) sprintf ( , _%03u'/ scene_num (i) ) ' ; ' ... 
deblank ( dataset (j,:) ) sprintf ( , _%03u', scene_num(j) 
num2str{ 100 * sidelap (count) ) ’%’ ]; 

if sidelap (count) > 0 

fprintf( fid, ' %s\n', line ); 
else 

fprintf( fid, 1 %s 1!! GAP !!!\n f , line ); 
end 
end 
else 

% flight lines in the opposite directions 
% estimate apparent endlap 


endlap_ = 1 - radius * ... 

abs ( ( vx(i) - vx(j) ) * ( longitude (j) - longitude (i) ) * ... 

cos ( ( latitude (i) + latitude (j) ) / 2 ) + ... 


( vy(i) - vy(j) ) * ( latitude (j) - latitude (i) ) ) / ... 
abs ( ( vx (i) - vx ( j ) ) * ( alongtrack (i) * vx(i) - alongtrack (j ) * vx(j) ) 


* 


* 


+ 


( vy(i) - vy(j) ) * ( alongtrack (i) * vy(i) - alongtrack (j ) * vy(j) ) ); 

if endlap_ > 0 

% include sidelap only when the endlap exists 

count = count +1; 

sidelap (count) = 1 - radius * ... 

abs ( ( vy(i) - vy(j) ) * ( longitude (j) - longitude (i) ) * ... 

cos ( ( latitude (i) + latitude (j) ) / 2 ) - 


vx( j) 
vy ( j ) 


( vx (i) - vx(j) ) * ( 
abs ( ( vx (i) - vx ( j ) ) * ( 


) + ... 

) ) ; 

line = [ ' >< 1 deblank ( 
deblank ( 
num2str( 100 * 
if sidelap (count) > 0 
fprintf ( fid, *%s\n', 
else 

fprintf ( fid, *%s !!! 
end 
end 


( vy(i) - vy(j) ) * ( 

dataset (i,:) ) sprintf ( 
dataset (j,:) ) sprintf ( 
sidelap (count) ) '%' ]; 

line ) ; 

GAP ! ! ! \n ' , line ) ; 


latitude (j) - latitude (i) ) ) / ... 
crosstrack (i) * vx(i) - crosstrack (j ) 

crosstrack (i) * vy(i) - crosstrack (j ) 

*_%03u', scene_num(i) )':'... 
'_%03u', scene_num(j) ) ' : ' ... 


* 


* 




s 

S 

l 

l 

i 
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end 

end 

end 

% Display final results 

fprintf ( fid, '\n%s\n\n', ident ); 

%fprintf( fid, 'Ground elevation: %g ft\n\n', elevation ); 


if exist ( 'endlap* ) == 1 


fprintf ( fid, 
fprintf ( fid, 
fprintf ( fid, 
fprintf ( fid, 
fprintf ( fid, 
fprintf ( fid, 

else 

fprintf ( fid, 
end 


' %d endlapsXn', length { endlap ) ); 

'Mean endlap : %f%%\n', 100 * mean ( endlap ) ); 

'Std. dev. : %f%%\n', 100 * std( endlap ) }; 

'Minimum : %f%%\n', 100 * min ( endlap ) ); 

'Maximum : %f%%\n', 100 * max ( endlap ) ); 

'Portion of endlaps smaller than %g%% : %f%%\n', ... 

100 * end_spec, 100 * length ( find( endlap < end_spec ) ) 

'No endlapsXn' ); 


fprintf ( 
if exist ( 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 

) ; 

else 

fprintf ( 
end 


fid, ' \n ' ) ; 

'sidelap' ) == 1 

( fid, '%d sidelapsXn', length ( sidelap ) ); 


( fid, 
( fid, 
( fid, 
( fid, 
( fid. 


'Mean sidelap 
'Std. dev. 
'Minimum 
' Maximum 


%f%%\n', 100 * mean ( sidelap ) ); 

%f%%\n', 100 * std( sidelap ) ); 

%f%%\n', 100 * min ( sidelap ) ); 

%f%%\n', 100 * max ( sidelap ) ); 


'Portion of sidelaps smaller than %g%% : %f%%\n', ... 

100 * side_spec, 100 * length ( find( sidelap < side_spec ) 


fid, 'No sidelapsXn' ); 


fclose ( fid ) ; 


/ length ( endlap ) ) ; 


) / length ( sidelap ) 
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